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Abstract 

For the two dimensional classical XY model we present extensive high - 
temperature -phase bulk data extracted based on a novel finite size scaling 
(FSS) Monte Carlo technique, along with FSS data near criticality. Our data 
verify that r/ = 1/4 sets in near criticality, and clarify the nature of correc¬ 
tion to the leading scaling behavior. However, the result of standard FSS 
analysis near criticality is inconsistent with other predictions of Kosterlitz’s 
renormalization group approach. 
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The classical XY model is defined by the Hamiltonian 

H = Y. ( 1 ) 

<ij> 

where s* is the 0(2) spin at site i and the sum is over nearest-neighbor spins. As a novel 
statistical model that undergoes a phase transition without long range order, the two di¬ 
mensional (2D) XY model is interesting in its own right, as well as being a model for 2D 
layers of either superconducting materials or films of liquid helium. 

The mechanism of this kind of phase transition was first illustrated by Berezinskii 
HI], and by Kosterlitz and Thouless [@](BKT) based on vortex binding scenario: At low- 
temperatures spin waves are the only significant fluctuations, and vortices are bound in 
pairs of zero vortices, thus not affecting the spin-wave description qualitatively; at high 
temperatures, however, the binding of the vortices decreases. The prediction 0 obtained 
with a renormalization-group method is that in the high temperature phase the critical 
properties of correlation length (,^) and magnetic susceptibility (y) are given by 


^(t) ~ exp(6t ^), 

(2) 

x{t) - 

(3) 


Here t = (T — Tc)/Tc, and the predicted values of z/ and rj are 1/2 and 1/4, respectively. 
For T < Tc, both ^ and y diverge identically with such a temperature dependent rj that 
r]{T) = 0. 

Kosterlitz’s renormalization group equations approach was extended to yield the correc¬ 
tions to Eq.(HD: 

{(()~exp(6r‘'^) ll + 0(()], (4) 

Accordingly x(t)/^(t)^/^ increases as t > 0. 

On the other hand, rigorous studies have established only two things: the existence of 
a phase transition at a finite temperature and also the existence of an upper bound of 
the critical temperature P|. Consequently the validity of Eq.(|D has been questioned both 
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analytically [0 and nnmerically p]-^. A different mechanism, namely “polymerization of 
the domain walls”, was proposed 0. Extensive Monte Carlo (MC) stndies [pA|-[T3[] np to 
^ ~ 70(1) and series expansions stndies snpported Eq.(2) over a pure power-law 

singnlarity. Bnt neither the MC data nor the series expansions showed r] = 0.25 and u = 0.5 
conclnsively. 

Primary difficnlties of various numerical methods in locating the parameters in Eq.(^ 
are due to four independent parameters involved so that, for example, a determination of 
u is sharply sensitive to the location of Tc, and vice versa. Thus, most MC studies have 


determined Tc by assuming z/ = 1/2 in Eq.(|), yielding 0.890 < Tc < 0.894 |1llHl^,PA|-^ 


The generic feature that has emerged from those MC studies was that the values of r] 
calculated from the thermodynamic values (TV value in the thermodynamic limit) of ^ and 
X are signihcantly larger than the predicted value |T^JT^JT6[| . 

In this work we present extensive MC data of the thermodynamic values of ^ and y up 
to T=0.93 (.^ ~ 1391(22)), indicating for the hrst time that rj = 1/4 holds as T —>• Tc. The 
unconstrained £t of our thermodynamic data are remarkably consistent with the prediction 
of Kosterlitz, yielding u ~ 0.48 and Tc — 0.893. The nature of the corrections to the leading 
scaling behavior Eq.(^ is clarified. We also attempt to locate Tc by conventional ESS , i.e., 
by hnding the temperature where xl holds with 1 ] = 1/A and, at the same time, 

where the 4th order cumulant ratio is invariant with L. This procedure yields Tc over 
the range 0.900 < Tc < 0.905, being incompatible with that estimated by assuming Eq.(2). 

Our extractions of TVs are based on a novel ESS technique |1T6tp2| -|2^ that facilitates 
drastically the MC measurements of TV. Namely, thermodynamic data can be computed 
from MC measurements on a much smaller lattice than required for the traditional direct 


measurements. The technique is based on the fundamental formula of ESS [If: 


Al (t) = Aoo (t) Qa (x), X = {t)lL. 


(5) 


Here Ai denotes a multiplicative renormalizable quantity A, defined on a hnite lattice of 
linear size L, and Qa is a universal scaling function of a different scaling variable ^L{t)/L. 
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Eq.(^ is supposed to hold even for the 2D XY model For a detailed description of 

the technique see Refs. |T^J2^J2^ , where the TVs extracted by employing the technique are 
shown to agree completely with those measured traditionally for the 2D and 3D Ising models 
23| , 2D XY model [0, and 2D Heisenberg ferromagnet p^ , |2^ . 

Standard theory of FSS predicts the presence of Tmm below which certain corrections to 
Eq. @ become non-negligible. It is shown for the 2D and 3D Ising models [^] that at least 
for L > 16 Eq.(D holds (hence, Lmin < 16) within the accuracy of the relative statistical 
errors of 0.3 percent. For T not so close to the criticality in the high temperature phase of 
the XY model, Eq.(||) holds with Lmm — 30 |^. As will be seen, however, Lmm increases 
as T —> Tc. Whether or not Tmm will increase endlessly as T —>■ Tc (thus violating FSS) is 
not clear at present, although its violation seems unlikely on theoretical ground ||^. We 
will assume that Eq.(^ is accurate for any value of t provided L > 80 . 


Here the extractions of the TVs of ^ and x are made in the range 0.93 <T < 0.98, based 
on available data of {x, Qa{x)) for T = 1.0 0. The values of L over the range 80 < L < 480 
are used for our analysis, and most interpolations are made in the range Qa{x) > 0.1. We 
employed Wolff’s single- cluster algorithm for the MC simulations, imposing periodic 
boundary condition on a square lattice. For a given set of T and L, typically 20 - 40 numbers 
of bins, (5 - 10 bins for the conventional FSS studies near criticality for the location of Tc), 
have been obtained for the measurements. Each bin is composed of 10 000 measurements, 
each of which is separated by 3-12 consecutive one cluster updating. Accordingly, the relative 
statistical errors of most raw-data calculated through the jack knife method are typically 
less than 0.1 percent for the x and 0.2 percent for the Correlation lengths are measured 
by the low-momentum behavior of the propagator 


22-24|. A comparison of our 


thermodynamic data in Table(l) with those obtained from a strong coupling analysis can 
be found in Ref. |^, which shows good agreement. 

Based on the bulk data, Qa{x) for A = .^ and x are plotted in Fig.(l), displaying an 
excellent data-collapse for L > 80. The data-collapse verihes Eq.(^) for the given values of 
L and T. The hgure also indicates a deviation from the data -collapse for the data point 
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corresponding to T = 0.95 and L = 60. Similar figures including data points with L < 60 
will be presented in a detailed paper |^. They will serve for a clear demonstration for the 
effect of correction to FSS for these values of L. Fig.(2) depicts xiT)/^{Ty~^ with r/ = 1/4 
as a function of ln(,^). It shows that the effective value of rj gradually decreases to r/ = 1/4 
as T ^ Tc- Empirical formula obtained from the data in Fig. (2), assuming = 0.893 (see 
below), is as follows: 


= <1+ *l ln(«)r, 


( 6 ) 


where a ~ 1.833, b = 0.955, and r ~ —0.413. Note that the slow decrease of x{i)/ 
as f 0 is inconsistent with any positive value of r such as was predicted in [Q. Without 
taking account of the corrections to Eq.(^, it turns out, however, that rj ~ 0.272(8). The 
value of x^ per degree of freedom {xYI^ df) for the latter case turns out to be approximately 
5.4 X 10^, which is much larger than that for Eq.(^. Accordingly, correction to Eq.(|^) is 
essential. 

We also employ the conventional FSS technique to locate independently. Namely, we 
measured both xl and the 4th order cumulant ratio defined as Ul = ?>—< >‘^ 

with S‘^ = Such well-known standard FSS properties at criticality as the scale 

invariance of the Ul and xl ~ are direct consequences of the fundamental FSS formula, 
Eq.(^. Therefore those properties should be valid only for L > L^in- We would like to 
stress, moreover, that an additional criterion should be satisfied in order for the former 
property to be true: a weak hyperscaling relation that the renormalized four point coupling 
at zero momentum remains a constant in the scaling region |^. Such a scaling 

behavior of has recently been shown numerically ||^ for the 2D XY and 0(3) vector 
models. 

We measured Ul and xl by varying L from 20 to 600 for some temperatures over the 


range 0.89 <T< 0.92. Certain non-asymptotic FSS behavior was observed for L < 80 [28 


which seems to be a consequence of the correction to Eq.(^ [^. Our data with L > 80 fit 
remarkably well to a standard FSS formula at the criticality 


5 






















( 7 ) 


XL ~ 

for T < 0.905 {x^/Ndp < 0.5) (Fig.(3)). The extracted values of rj from the linear fit 
are: 7]{T) = 0.260(1), 0.252(0), 0.245(0), and 0.231(1) for T=0.91, 0.905, 0.90, and 0.89 
respectively. Thus r]{Tc) = 0.25 is consistent with 0.90 < Tc < 0.905, most probably 
Tc ~ 0.904. This estimate is also consistent with our data of Ul (Fig.(4)). Ul{T = 0.92) 
definitely decreases with L from L=80 to L=240, so that Tc is obviously smaller than 0.92. 
At T=0.91 and 0.905 it varies so mildly that some fictitious scale invariance of Up appears 
over a range of L not too large; nevertheless, it eventually decreases. This fictitious scale 
invariance may result from an extremely large value of ^ over the range of temperature. We 
observe the invariance of UpiT) with respect to L at T = 0.90, indicating Tc > 0.90. 

We carried out a fit of ^ data over the range 5.01(3) < ^ < 1391(22) (corresponding 
to 0.93 < T < 1.19) to the exponential singularity, Eq.(||). Our fit, based on down-hill 
simplex algorithm, is highly non-linear for Eq.(^), so that we have treated Tc as an input 
parameter over a reasonable range. The values of x^/Ndf and u as a. function of the 
input Tc are plotted in Fig.(5). The best £t is obtained for Tc ~ 0.893 and for u ~ 0.48 
{x^/Ndf ^ 0.75). u = 0.5 holds at Tc ~ 0.892 with x^/Ndf — 0.80, being in agreement with 
the range of Tc conventionally accepted. The agreement is gratifying and may be regarded 
as another verification for Eq.(|^) with = 80, at least for T > 0.93. The results of the 
unconstrained fit agree extremely well with the prediction in Eq.(2). Notice, however, that 
T] = 1/4: holds for Tc — 0.904 according to the conventional FSS, being much larger than 
0.892. Assuming Tc = 0.904 the best fit is for f ~ 0.30(1) with x^/Ndf — 6.69; assuming 
both Tc = 0.904 and z/ = 0.5 the best fit has x^/^df — 263. Therefore, u = 1/2 in Eq.(2) 
is definitely inconsistent with the prediction of FSS at criticality. 

We also checked whether or not the data fit to ^ = ci(l -t- C 2 t^) exp(6t“‘') with the input 
value of Tc over 0.90 < Tc < 0.905. It turns out that they fit to the formula for some 
negative value of 6 only (being inconsistent with Eq.(4)), with the value of x^/Ndf being 
almost insensitive to the choice of Tc. However, this function has too much freedom to the 
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fit and gives unstable predictions for the values of both Tc and u. 

On the other hand, the data £t very well to a modihed second order phase transition 
(see also Ref. H). 

^ = Ci{l + C2t^)t ( 8 ) 


Assuming Tc = 0.904 we estimate that z/ ~ 3.10(8), 9 ~ 1.82(8), and C 2 — 44.2 with 
X^/Ndf — 1-2 (Fig.(6)). As the value of the hxed Tc becomes larger, the data £t better 
to the modihed power-law singularity. The large value of C 2 is unusual, implying that at 
sufficiently large t the correction term becomes dominant with an effective exponent of ^ 
being equal to u — 6. We note that our effective u for large t agrees with the prediction in 
E3- For t small enough the correction becomes negligible regardless of the value of C 2 , while 


over the intermediate t the effective critical exponent changes gradually. 

Assuming Tc ~ 0.904, even a modihed exponential singularity including correction is 
unlikely: For a power-law singularity, ^ ~ t~’^, the FSS of Uiit) is given by, UL{t) ~ 
fu{L^^^t) [^. Using the analyticity of Ul for a hnite L, after a Taylor expansion we 
get Ui(t) ~ Ul{0) -|- for a sufficiently small value of the An 

exponential singularity corresponds to z/ —cx) so that Ul(T) at T = 0.905 (t ~ 10“^), 
for example, would be L independent regardless of the value of L, contrary to our hnding 
in Fig.(4). Similarly, if an exponential singularity is valid Xl{T) l2-ri(T,) Yl0\d 

identically at any T close to T^; in other words, rj would remain unchanged over a hnite 
range of T close to T^, which again contradicts our data in Fig. (3). In view of this FSS 
argument, an exponential singularity could be consistent with our data of xl (Fig-(3)) and 
Ul (Fig.(4)) only if Tc is considerably smaller than 0.904. Assuming a power-law singularity, 
on the other hand, one can easily check from the expansion of Xiit) that even for such a 
modestly large value of z/ as z/ ~ 3.1, ln(xi/LF4^ T=0.905 varies very mildly with L in 
agreement with our data in Fig.(2) p8| . 

It remains to be seen, however, whether or not Tmm diverges as T —T^. In this case, 
the conventional FSS analysis at the criticality would be misleading. It is also possible that 
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the correction to the leading scaling behavior (Eq.(|^)) resnlts in a modihcation to the FSS 
behavior (at criticality). The conventional analysis then would not give a very accurate 
estimate of r}{Tc). A more detailed account for the possible failure of the standard FSS 
analysis at the criticality, along with details of our analysis and methods of obtaining our 
data, will be presented in a longer paper . 


In conclusion we have reported strong numerical evidence that rj = 1/A sets in at temper¬ 
atures close to the criticality. The nature of the correction to the leading scaling behavior 
Fq.(^ has been clarihed. The unconstrained £t of our high temperature thermodynamic 
data yields results consistent with the predictions of the BKT for the values of u and r]. 
This, nevertheless, does not rule out the possibility of the ordinary second order phase 
transition. Knowing the true critical point is crucial to the resolution of the order of the 
phase transition in the 2D classical XY model. The critical point obtained from the conven¬ 
tional FSS analysis, however, seems to be incompatible with the conventional exponential 
singularity. 

After the completion of the current work we have become aware of a recent related work 
by Kenna and Irving that questions the vortex binding scenario in the 2D XY model 


based on FSS study of Lee-Yang zeros. They argue that Fq.(2) cannot be compatible with 
T] = 1/4: without a certain correction. However, the proposed correction, r ~ 0.02 in Eq.(^), 
is inconsistent with, our picture in Fig.(2) in that their x{t)/increases as t ^ 0. 

The author would like to thank Prof. M. Howard Lee for his proof reading of this 
manuscript. Current work was partially supported by NSF grant DMR-9405018. 






TABLES 


TABLE I. The thermodynamic values of ^ and x for 0.93 < T < 0.98, extracted through our 
MC ESS technique. 

T 0.98 0.97 0.96 0.95 0.945 0.94 0.935 0.93 

C 70.4(4) 100.3(7) 155.7(1.5) 262.7(2.9) 364.5(3.6) 539.3(4.5) 846.7(6.4) 1391(22) 

X 4284(20) 7932(30) 16978(55) 42295(108) 74742(206) 147536(614) 3.24(1)xlO^ 7.70(9)xlO® 
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Figure Captions: 

Fig.(l): Q^{x) (upper data) and Qxi^) (lower data) for L > 60, calculated with the data 
from Table(l). It should be clear that for each T a data point with a larger L corresponds 
to a smaller x. For T = 0.95 and L=60, the deviation from data-collapse is apparent. For 
the sake of the higher resolution of the data, we show only a lower part of Q^{x)-, for the 
upper part, see Ref. []^ . 

Fig.(2): versus ln(,^(T)) for T over the range, 0.93 <T < 1.25. The thermo¬ 

dynamic data for T > 0.98 are taken from Ref. |^, which agree completely with the data 
in [^ [Q. The tendency to decrease implies rj > 0.25, but the values of tend to 

stablize at approximately In(^) > 5. 

Fig.(3): ln(x/,/L^/^) versus ln(L) over 80 < L < 600. Here the slope of a straight line 
is related to the value of rj, that is, g = 1/4 — slope. Note a slow decrease (increase) of 
ln(xi/L’^/^) with respect to ln(L) for T=0.905 (T=0.90); hence, 0.90 < < 0.905 with 


g = 1/4. 

Fig. (4): Ul as a function of ln(L) over 80 < L < 600, for various temperatures near 
criticality. Ul is L independent (within small statistical error) over 80 < L < 150, and 
80 < L < 360 for T = 0.91 and 0.905 respectively (pseudo scale invariance), although it 
eventually decreases with L. 

Fig.(5); The values of x^I ^df (lower) and v (upper) as a function of the input parameter 
Tc, which are obtained from the x^ fit of the ^ data to Eq.(^. 
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Fig.(6): In(^) as a function of | lnf| assuming Tc = 0.904. The circle symbols denote our data 
in Tab.(l), while the smooth curve represents the modihed second order phase transition, 
Eq.(^, with the values of the parameters reported in the text (except for ci ~ 2.19 x 10“^). 
The error bars are almost invisible. Here a straight line would represent a pure power-law 
singularity, with its slope equivalent to the value of u. Note that for | lnf| > 2.6 the data 
are almost linear, implying that the effect of correction is already not so important in this 
regime. 
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